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GLOBAL SOLUTIONS FOR A SUPERCRITICAL 
DRIFT-DIFFUSION EQUATION 


JAN BURCZAK AND RAFAEL GRANERO-BELINCHON 

Abstract. We study the global existence of solutions to a one-dimensional 
drift-diffusion equation with logistic term, generalizing the classical parabolic- 
elliptic Keller-Segel aggregation equation arising in mathematical biol¬ 
ogy. In particular, we prove that there exists a global weak solution, if 
the order of the fractional diffusion a € (1 — ci,2], where ci > 0 is an 
explicit constant depending on the physical parameters present in the 
problem (chemosensitivity and strength of logistic damping). Further¬ 
more, in the range 1 — C 2 < Of < 2 with 0 < C 2 < ci, the solution is 
globally smooth. Let us emphasize that when a < 1, the diffusion is in 
the supercritical regime. 
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1. Introduction 
T he drift-diffusion equation 

(1) dtu =—vk°‘u+ 'S/■ {uB{u)) + f{u), 

where B{u) is typically a vector of nonlocal operators and A = y/—k. (see 
(11) below), appears widely in applications. The parameter 0 < a < 2 is 
the order of the diffusion and it measures the strength of the viscous effects. 

For instance, the two dimensional incompressible Navier-Stokes equations 
in its vorticity formulation can be written as 

( 2 ) dtu = ku+ V ■ {uV^ 

where V"*- = {—dx 2 )dx-^)- This equation governs the motion of two-dimensional, 
incompressible, homogeneous fluids in absence of forcing (see [48]). Equa¬ 
tion (2) can be recovered from equation (1) by taking u = l, a = 2, / = 0 
and B = 

Another famous equation akin to (1) is the parabolic-elliptic simplification 
of the Keller-Segel system with logistic source 

(3) dtu = An -|- V • (nVA~^n) + — ru^, 

It appears as a model of chemotaxis, i.e. the proliferation and motion of 
cells (see the pioneer work of E. Keller & L. Segel [41] and the reviews by 
A. Blanchet [9] and Hillen & Painter [40], whose cell-kinetics model (M8) is 
the doubly-parabolic version of (3)). Here n > 0 is the density of cells. To 
obtain (3) form (1) one takes u = 1, a = 2, B = VA“^ and f = fjLU — rv?. 
Equation (3) is also a model of gravitational collapse (see the works by Biler 
[5] and Ascasibar, Granero-Belinchon &: Moreno [1]) 

Eurthermore, in one spatial dimension, the equation 

(4) dtu = -dx{uHu), 

where H denotes the Hilbert transform (see (10) and (12) below), has been 
proposed as a model of the dynamics of a dislocation density u (see [31] and 
the work by Biler, Karch &: Monneau [7]). Equation (4) appears also in a 
totally different context, namely as a one-dimensional model of the surface 
quasi-geostrophic equation (see Castro &: Cordoba [16]). In order to recover 

(4) from (1), we choose u = 0, / = 0 and B = —H. 

Einally, notice that the famous Burgers equation 

(5) dtu = —A"u — udxU 

can be obtained from (1) by taking B(u) = —uj2. 

In all these equations there is a competition between the diffusion term 
given by 

A“ 

and the transport term 

V • {uB{u)). 

Consequently, if the order of the diffusion a is large enough, the solution u 
will remain smooth for every time. On the other hand, when the order of 
the diffusion is too weak, there are finite-time singularities. 

Burgers equation (5) is the paradigm here. There, if a > 1 (this is known 
as the subcritical case), there exists a global, smooth solution, whereas for 
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a < 1 (this is known as the supercritical case), there is a finite time T*, 
such that 

limsup \\dxu{t)\\L^ = oo. 
t->T* 

The remaining case a = 1 is known as the critical case. Here the diffusive 
operator is of order one, so it perfectly matches the order of the transport 
term. Kiselev, Nazarov &: Shterenberg [43], Dong, Du & Li [33] and Con¬ 
stantin &: Vicol [22] proved that in the critical case, there exists a global, 
smooth solution. Furthermore, Dabkowski, Kiselev, Silvestre & Vicol proved 
the global existence of classical solutions for a logarithmically supercritical 
Burgers equation [30]. 

A similar picture appears for the one-dimensional case of (1) with / = 0, 
which is a fractional-diffusion version of Keller-Segel system (3) with // = 
r = 0. This is because the Burgers equation (5) is a primitive equation to 
this one. For more details on available results, see Section 2. 


1.1. Considered problem. In this paper we are going to study the drift- 
diffusion equation (1) on the one-dimensional torus T, i.e. 

( 6 ) dtu =—A^u + xdx{uB{u)) + f{u), xGT 
with 

(7) B{u) = Af^-^H{l + A^)-\, 
and / of logistic form 

(8) f{u) = ru{l — u), 0 < r. 

Our primary motivation here is the (parabolic-elliptic) Keller-Segel model, 
which corresponds to (3 = 2. 


1.2. Outline of results. Our results concern the global in time existence 
of solutions to (6) with a < 2 and x > f. Primarily, we will show that for 
positive 

a > 1 - r/x, 

equation (6) possesses a global smooth solution, see Theorem 2. To the 
best of our knowledge, this is the first result on global smooth solutions in 
the supercritical regime of an equation modelling chemotaxis. This result 
may be somehow surprising, because the global existence is obtained due 
to an inviscid regularization —rv?, present in the logistic term that may be 
arbitrarily small, i.e. 0 < r <C 1. Let us recall that for r = 0 the critical 
case a = 1 is globally well-posed, see the work [12] by the authors. 

Along with the smoothness result, we also study the global existence of 
weak solutions and show that there exists a global weak solution for 


a > 



for x/2 > r 
for x/2 < r 


see Theorem 1. At the cost of reducing the integrability and smoothness of 
our weak solutions, we can relax the condition on a and obtain 


a > 


) 1 —s — 
' 2-l-s 


for x/2 > T and 0 < s < 
for x/2 < r. 
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This result is more difficult to obtain than its analogue for the Burgers 
equation (5), due to lack of the energy balance. These main theorems 
cover the case of an arbitrary, finite time of existence T (with no small¬ 
ness assumptions on data). They are supplemented by Proposition 1 on 
boundedness of u for T = oo. It holds for supercritical or critical diffusions, 
where in the latter case we need an explicit smallness condition on data or 
chemosensitivity parameter y. 

Let us remark here that our result may be easily adapted for another class 
of operators B. For instance, one may consider the case 

(9) B{u) = - {u)), 


where 


(u) = — [ u{x)dx. 
2vr Jf 


Moreover, the same ideas work for the case of the spatial domain being 


1.3. Plan of the paper. In Section 2 we review some results on drift- 

diffusion equation that are related to our problem. Next, in Section 3, we 
state our main results. In Section 4 we provide some preliminary results, 
including the local existence and analyticity and certain estimates for lower 
order norms. Section 5 is devoted to the proof of Theorem 1 on global in 
time existence of weak solution, whereas in Section 6 we prove Theorem 
2 on global in time existence of smooth solution. In Section 7 we prove 
Proposition 1 on boundedness for T = oo. Finally, in the appendix we 
provide some estimates for the fractional laplacian and inequalities for the 
standard Llog(L) entropy and the seminorm. 

1.4. Notation. We write H for the Hilbert transform and A = \/—A, i.e. 

( 10 ) Hu{^) = 

(11) A^(o = ierfi(o, 

where " denotes the usual Fourier transform. Notice that in one dimen¬ 
sion A = dxH and Hu{0) = 0. These operators have the following kernel 
representation 

(12) «„(,) = ip, 


(13) 

We also have 

(14) 
with 


Au{x) = —P.V. 
27r 


u{x) - u{y) 


't sm 




K^u{x) ^ Cc,( ^ f 


u{x) — u{x — ri)dri 
u{x) — u{x — rj)drj 


Cry - 


+ P V . 

Jt 

F(1 -|- a) cos((l — Q;)7r/2) 

TT 


> 0 . 


We denote by TA the periodic heat kernel at time t = e. 
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We write T^ax for the maximum time of existence of a given solntion u. 
Finally, to simplify the notation, for the case B{u) given by (7) we define 

(15) n = (1 + 
while in the case (9) 

(16) V = A~^{u — (u)). 

Notice that, in both cases, v can be defined using Fourier series. Further¬ 
more, we have 

B{u) = 


1.5. Functional spaces. The fractional L^-based Sobolev spaces, W^’P{T), 


are 


W^’P = (/ G LP(T),4"V G LP{T), e LP{T X T)1 , 


\x - y[- 


+(«-LsJ) 


with norm 


p _ 


LP 


+ 


p 

W^’P' 


p _ 


VIIlp + 


T JT 


\dxhix) - dlf^f{y)\P 

\x — y|l+(«-W)P 


dxdy. 


In the case p = 2, we write F7^(T) = ly^’^(T) for the standard non- 
homogeneous Sobolev spaces with norm 

\ h ^ = IIA^IIl- 


2 , II rii2 

L 2 + 




2. Prior results for one-dimensional drift-diffusion equations 

In this section we are going to review some results concerning certain 
one-dimensional models. 

First, let us recall the following nonlocal version of the Kuramoto-Shivashinsky 
equation: 

(17) dtu + udxU + eA^~^^u = A'^u, 0 < 7 < 1-|-<5, e,(5>0. 

It provides a simple model for the stabilization of a Hadamard instability, 
with growth rate proportional to the absolute value of the wavenumber, 
by higher-order viscous diffusion. The Kelvin-Helmholtz instability for the 
Euler equations is an instability of this type. Equation (17) was studied 
by Granero-Belinchon & Hunter [37]. These authors showed the chaotic 
character of the dynamics and studied analytically some properties of the 
attractor and the solutions. 

Next, let us recall that extensive studies were performed for the so called 
(5—model 

(18) dtu + (1 — 5)HudxU + 5dx{uHu) + uA"u = 0, 0 < (5 < 1. 

Let us for a moment focus on the the inviscid case u = 0. This equation 
appears as a one-dimensional model of the inviscid surface quasi-geostrophic 
(SQG) equation. From the viewpoint of applications, the SQG eqnation 
models the dynamics of hot and cold air in the atmosphere; from a purely 
mathematical perspective, it is widely studied as a two dimensional model 
of the three dimensional Euler equations in its vorticity formulation. Let 
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US note that there are several singularity formation results available. The 
singularity in the cases 0 < 5 < 1/3 and <5 = 1 was proved by Morlet [49]. 
In the range 0 < 5 < 1 the finite time singularity was proved by Chae, 
Cordoba, Cordoba &: Fontelos [19]. Finally, the singularity in the case 5 = 0 
was proved by Cordoba, Cordoba and Fontelos [26] and by Silvestre &: Vicol 
[52]. The case 5 = — 1 and a > 1 is similar to the nonlocal Kuramoto- 
Sivashinsky equation (17) and it has been studied by Li &: Rodrigo in [47]. 
On the other hand, Bae &: Granero-Belinchon [3] proved the global existence 
of weak solution in the case 1 > 5 > 0.5. 

The two extreme cases 5 = 0 and 5 = 1, both in the inviscid and viscous 
case, deserve special further attention. 

The case 5 = 0 reads 

(19) dtu + HudxU + uA^u = 0. 

This is a Burgers-like equation with a nonlocal velocity. It is contained in 
the scale of equations 

(20) dfU + A^ HudxU + i^A°‘u = 0, sE(—1,1). 

The family of equations (20) was proposed as a model of the two-dimensional 
a—patches by Dong &: Li [34]. 

As we have outlined, in the inviscid case, equation (19) develops finite 
time singularities. On the other hand, one can prove that if a > 1 there 
exist global classical solutions. The global existence for arbitrary initial 
data in the critical case a = 1 was proved by Dong [32], see also Kiselev 
[42]. The singularity formation in the viscous case has been proved in the 
range 0 < a < 0.25 by Li &: Rodrigo [46]. The range 0.25 < a < 1 remains 
a major open problem. In this regard, it is worth mentioning that Silvestre 
&: Vicol [52] conjecture that solutions of (19) with 0.25 < a < 1 which arise 
as limits from vanishing viscosity approximations are bounded in the Holder 
class 1/2. 

The case 5 = 1 of (18), or, equivalently, equation (4) in the inviscid case, 
has been studied by many authors. For instance, global classical solutions 
for initial data satisfying a strict positivity condition and finite time singn- 
larities otherwise are proved by Castro and Cordoba in [16] (see also [17]). 
The existence of global weak solutions is proved by Carrillo, Ferreira &: Pre- 
cioso [15] or Bae & Cranero-Belinchon [3]. The authors studied in [11] the 
asymptotic behavior. In particular, we obtained the following inequality 

— fv(rt)--X(rt)^ > 0, 

dt \ 2mm2;Uo(x) J 

where the entropy and the Fisher information are dehned as 

F{u)= / u\og{u) — udx, 

JT 

I{u) = f dx, 

Jr 

respectively, and uq such that (uq) = 1 denotes the initial data for equation 
(4). As a consequence, we have proved the following nonlocal logarithmic 
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Sobolev inequality 

( 21 ) f fix) logifix))dx < 277 + ^ . f \A^-^fix)\^ dx. 

JT ^ miiix I { X ) Jj 

Inequality (21) remains valid for every function f G n H^'^, / > 0. 
Consequently, for strictly positive and bounded below data we have obtained 
exponential decay of J^(u(f)). 

Equation (4) is a particular case of the aggregation equation 

(22) dtu = dxi-fiiu)Hu + xudxidl)~^iu - (u)) + rtt(l - u), 

studied by Granero-Belinchon Sz Orive-Illera [39] and the authors in [11], In 
turn, equation ( 22 ) is akin to the classical (parabolic-elliptic) Keller-Segel 
equation in one spatial dimension: 

(23) dtu = -A°‘u + dxiud^{dl)~^{u - (u))). 

Since Keller-Segel systems are one of our main interests, let us concentrate on 
it for a moment. Equation (23) is the one dimensional model of the behavior 
of microorganisms (with density u), where their motility follows A", instead 
the (more classical) —d^. This choice for the diffusive operator is supported 
by a strong evidence, both theoretical and experimental, that feeding strate¬ 
gies based on Levy process are used by organisms in certain situations. For 
instance, the interested reader can refer to Lewandowsky, White &: Schus¬ 
ter [45] for amoebas, Klafter, Lewandowsky & White [44] and Bartumeus, 
Peters, Pueyo, Marrase &z Catalan [4] for microzooplancton, Shlesinger & 
Klafter [51] for flying ants, Cole [20] for fruit flies, Atkinson, Rhodes, Mac¬ 
Donald & Anderson [2] for jackals and Raichlen, Wood, Gordon, Mabulla, 
Marlowe, &: Pontzer for the Hadza tribe [50]. 

Equation (23) was hrst studied by Escudero [35]. He proved the global 
existence of solution in the case 1 < a < 2. This result was later improved by 
Bournaveas & Calvez [10], where the authors proved hnite time singularities 
for the case 0 < a < 1 and the existence of AT > 0 such that for the case 
a = 1 and the initial data satisfying the smallness restriction 

||woIIli < 

there exists a global smooth solution. Ascasibar, Granero-Belinchon & 
Moreno [I] obtained the global existence in the case a > 1 and a = 1 
under the relaxed smallness condition restriction 

ll-olii. < T. 

Under this restriction, the solution u{t) verifies 

lln(t) - (rio)||L°° 0 . 

All these results where improved by the authors in [12], where we showed 
the global existence and smoothness in the case a = 1 and arbitrarily large 
initial data. Furthermore, we also showed that if the initial data uq is such 
that 

IImoIIli < 277, 

then the solution tt(f) verihes the asymptotic behaviour 
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for certain S, cr > 0 . 

Let us mention also here the works by Biler & Karch [ 6 ], Biler &: Wu [ 8 ] 
and Wu &: Zheng [53], where nonlocal Keller-Segel-type models with very 
similar in spirit operators B{u) are considered. 

Finally, let us mention in our review of one-dimensional models the equa¬ 
tion 

(24) = -A« + 8, { 0 ^) ■ 

It arises as a model of the interface in two-phase, porous flow. Equation 
(24) was considered in [28] by Cordoba, Gancedo & Orive-Illera, where the 
authors studied the local existence and qualitative behavior of the solutions. 
Granero-Belinchon, Navarro &: Ortega [38] proved the global existence of 
weak solutions. 


3. Statement of results 


We start with our definition of weak solution 


Definition 1. u{x,t) E L^(0, T;L^) is a global weak solution of ( 6 ), (7) (or 
of ( 6 ), (9)j with logistic forcing ( 8 ) if for allT > 0 and (f E P([—1,T) x T) 
we have 



u{—dt4> -|- A“(/)) -|- uB{u)dx'f — f{u)(l)dxds 


uo4>{0)dx. 


Let us observe that the bigger the ratio y/r is in ( 6 ), the stronger is 
the destabilization coming from the term xdx{uB{u)), compared with the 
logistic damping f{u) = ru{l — u). Nevertheless, we will be able to show 
that for any ratio x/r, there is certain region in the supercritical regime, 
where the stabilizing effect prevails. 

Our first result is the global existence of weak solutions with certain higher 
regularity. Define 

/ ^ 

s = min -, 1 

Vx - J' 

and s~ as any number strictly smaller than s. 


Theorem 1. Take a, x,/3,r > 0. Let uq > 0, uq E be the initial data, 
(i) Case r < x/2. 

Assume that 


(25) a > 1 - s. 

Then, there exists a global in time weak solution 

u E l2+"“(0,T;L2+"”) 

for ( 6 ),(7) (or ( 6 ), (9) j with ( 8 ) (in the sense of Definition 1). This solution 
verifies for any T < oo 

u E L°°(0, T; L^+*) nL2(0, T; T; pp^«/(2+2s)-5,i+8)^ 

where 0 < 6 < 0^ and 0 < di < ^ are arbitrarily small numbers. 

At the cost of decreasing s to s~ above, one can relax condition (25) and 
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assume instead 

1 _ s _ q 2 

( 26 ) 

(a) Case r > x/2. 

In addition to the previous case, we have here for any T < oo 

In other words, for arbitrarily large ratio x/r we have global existence of 
weak solutions in the supercritical regime a > 1 — s. Moreover, once the 
logistic term is strong enough (case r > x/ 2 ) every diffusion a > 0 is enough 
to provide us with global solutions. This result is particularly interesting 
for the case 0 < r ^ 1 . 

Remark 1. Notice that in the case r = 0 and every a > 0, global existence 
of fat least one) weak solution in the measure space H~^ can he obtained 
(see [12]). 

Our second result establishes the global existence of smooth solution for 
a range of diffusions max(l — r/x, 0 )<a <2 and sufficiently smooth initial 
data. In particular, notice that for any r > 0 there are diffusions with 
order less than one (supercritical range) such that the solution is global and 
smooth. 

Theorem 2. Take r,x,o:, /3 > 0 and uq > 0, uq € Then, if 

r 

1-< a < 2, 

X 

there exists the unique global in time solution for ( 6 ),(7) (or for ( 6 ), (9) J 
with ( 8 ), verifying for any T < oo 

uGCi[0,T],H^). 

Furthermore, u{t) is real analytic for any t G (0, T] and verifies 
sup ||u(t)||^i < max{27r, ||uo||li}' 

0<t<oo 

In particular, we obtain the following corollary: 

Corollary 1. Given r > x > 0, the system (6), (7) (or (6), (9)) with (8) 
posesses a global in time, classical solution for every initial data uq G H^. 

The idea of the proof of Theorem 2 is to use pointwise bounds on the 
A“u term that control the reaction term and prevent blow-up. This method 
has been developed by Cordoba &: Cordoba [24, 25] and Constantin &: Vicol 
[23]. It has proved itself useful in many different equations: see [25, 22] for 
the critical surface quasi-geostrophic equation, [27, 28, 18, 29, 21, 36] for its 
application to the Muskat problem and [1, 39, 11] for its application to the 
aggregation equations of Keller-Segel type. 

Notice also that the solutions are global but, as far as we know, they may 
be growing in ||M(t)||j,<x) without bound. In other words, the logistic term is 
not strong enough to provide global, uniform estimates. In particular, the 
estimates for ||n||i,cx) depend on the estimate for the diffusion term. In this 
context, let us provide a result concerning the large-time behaviour: 
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Proposition 1 . Let r,x, f3 > 0, a > 1 and uq > 0, uq G H^. Then, the 
global in time solution for (6), (7) (or (6), (9)) with (8) verifies 

sup \\u{t)\\L-o < C{a,r,x,\Wo\\L--,\Wo\\L^)- 

0<t<oo 

For a = 1 we have 

sup ||M(t)||L°o < C{r,x,\Wo\\L^,\\uo\\L^), 

0<t<oo 

provided x<r+ 2 ^„,ax(||lolLi ,2^) ' 

Corollary 2. For 0 = 1 the global in time solution for (6), (7) (or (6), 
(9) ) with (8) verifies 

sup ||?i(t)||L°o < C{a,r,x,\\uo\\L<^,\Wo\\Li)- 

0<t<oo 

if one of the following conditions is satisfied 

• X < r and the initial data may be arbitrarily large. 

• < 2-71 and x < r + In particular, for x < any r > 0 

is admissible. 


4. Preliminary results 

4.1. Elliptic estimates for v. We need certain elliptic estimates of v, de¬ 
fined in (15) and (16): 

Lemma 1. For suffieiently smooth, nonnegative funetion u and v given by 
(15) one has 

(1) min 3 ;i;(x) > 0, in particular A^v{x) < u{x), 

(2) ||u(t)||Loc < ||u(t)||ioo, 

( 3 ) \\A^v{t)\\L^<4\\u{t)\\L^, 

(4) \\ da , Al ^ v { t )\\ L <^ < 4 :\\ dxu { t )\\ L <^, 

(5) l\\A^+^v{t)\\l, + \\Af^/^+Mt)\\h < Vs > 0. 

Proof. Recall that v given by (15) means that it is a solution to 

(27) V + A^v = u. 

Part (1) Let us evaluate this equation at xt such that miua; v{x, t) = v{xt, t). 
Due to the definition of xt, we have v{xt) — v{xt — ??) < 0, and, consequently, 
using (14), A^i;(xt) < 0. Using this fact we get 

(28) v{xt) > A^v{xt) + v{xt) > minn(x,t) > 0 . 

X 

The non negativity of u starting from nonnegative data follows from a similar 
technique of tracking minimum, compare [ 11 ]. 

Part (2) In the same way, if we evaluate at xt such that imiXxv{x,t) = 
v{xt,t), we have A^v{xt) > 0 and we get the bound 

v{xt) < v{xt) + A^v{xt) < maxrt(x,t). 

X 

Part (3) and (4). Now we evaluate the equation (27) at Xt such that 
TaaiCx A^v{x,t) = Al^v{Xt,t) > 0. Due to part ( 2 ), we have 

(29) Af^v{Xt) = u{Xt) - viXt) < 2 ||u(t)||Loo. 
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Similarly, for Xt such that miiia, A^v {x,t) = Al^v{Xt,t) < 0 we have 

(30) -A^viXt) = -u{Xt) + viXt) < 2||u(t)|Uoc. 

Collecting estimates (29) and (30), we obtain (3). In the same way, we can 
prove (4). 

Part (5). It follows from (27). □ 

In turn, the case where v is dehned by (16) is studied in the following 
lemma. 

Lemma 2. For sufficiently smooth, nonnegative function u and v given by 
(15) one has 

(1) ||A^u(t)||Loc < \\u{t)\\Loo, 

(2) \\dxA^v{t)\\L-° < \\dxu{t)\\L^, 

(3) l\\A^+^v{t)\\l,<l\\AMt)\\h, Vs>0. 

Proof. In this case, u is a solution to 

(31) A^v = u — (u). 

Now, the proof is a simpler version of respective part of proof of the previous 
lemma. □ 

4.2. Local existence. We start with the local existence result. 

Lemma 3. Let uq G he a non-negative initial data. Then, i/O < a < 2, 
there exists a time T*{uo) > 0 such that there exists a non-negative solution 

u{t)eC{[0,T*{uo)],H^) 

to the equation (6), (7) (or to (6), (9) j with (8). Moreover, if for a given T 
the solution verifies the following bound 

u{s)\\Lacds < oo, 

then the solution may be extended up to time T -\- 5 for small enough 0 < d. 
Furthermore, under the restriction 1 < a < 2, the solution becomes real 
analytic. 

Once we have Lemmas 1 and 2, the proof of this result is similar to the 
proofs in [1, 39, 11, 13]. For brevity, we omit the standard proof. 

4.3. Auxiliary energy estimates. In the next lemma we collect some 
bounds for Lebesgue and Sobolev norms of the solution. 

Lemma 4. Let uq € be a non-negative initial data. Take s G and 
X > r > 0 such that 

r 

0 < s < -. 

X-r 

Given an initial data uq for equation (6), (7) (or for (6), (9)) with (8), let 
us define 

(32) AA = max{ II uo 11^1,271}. 

Then, the non-negative solution u{x, t) verifies for all times 0 < t < T < 
Tmax the following bounds 

(1) ll^(^)lli,i 
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(2) j^\\u{s)\\\2ds <Mt + 2M, 

(3) \\u{t)\\Ls+i < e’'*||uo||L»+i, 

(4) (r(s + 1) - xs) Jo lk(s)|llt+ 2 ds < 

(5) /o <C(a,s,5J|no||L^^^ 

for arbitrary 0 < 5 < a /(2 + 2 s), 

( 6 ) Jo Ms)\\'^^^j2-s,ids < C{a,6,J\f, ||? io || li +'*, 7 ") 
for arbitrary 0 < 6 < a/2. 

Proof. Part (1) and (2). The proof of this result is analogous to the proof 
in [13, Lemmas 1, 2]. 

Part (3). Testing ( 6 ) with we get after integrations by parts 
(33) 

^--f- [ u^~^^dx + [ u^A°'udx= [ (u^~^^)A^v + ru^~^^ — ru^~^^dx. 

IdtJj A JtS + 1^ ^ 


s + 

Lemma 7 implies 


0 < / u^(x)A°‘u(x)dx. 

Jt 


This and A^v < u (valid due to definition (27) and Lemma 1 or due to 
definition (31)) gives 


(34) 


A 

dt 


,S+1 


^S + 1 _|_ + 1 ) — xs) [ < r{s + 1 ) / 

Jt Jt 

Therefore, for s such that (x — t)s < r we get 
(35) ||u(f)||£,.+i < e'’*||uo||L-+i- 

Part (4) Integrating (34) in time we obtain that for any t G [0,T], 

(r(s + l) —xs) [ [ u^^^(x, s)(ix(is < [r(s + + 1] f Uo'^^{x)dx, 

Jo Jt Jt 

where the latter comes from (35). 

Part (5) Integrating (33) in time and using (35) together with Lemma 
7, we obtain 

rT 


\u 


|2+2s 

I ^X“/(2+2s)-i5.1+s 


< C{a,s,6,uo,T). 


Part (6). We define the functional 

-F = / Mlog(M) — ti + Idx. 


We have 


—iF = [ 9frtlog(u)(ix 
dt Jt 


< - 


< - 


/ A"rilog(rt)dx + X / A^vudx 
Jt Jt 

/ A"ulog(u)(ix + xlh 
Jt 


iMh- 
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Consequently, we obtain 


Using 


J^{u{t)) + / / A'^ulog{u)dx < J^{u{0)) + Mt + 2J\f. 

Jo Jt 

T{u{0)) < / tt(x, + 2-71 < ||n(0)||£,i+s + 271, 

J {ti(ai,0)>l} 

(36) T{u{t)) + f [ A“(n + 1 ) log(tt + 1 ) < ||uo||li+'> + 2 vr + 2 xAA(r + 2 ). 
Jo Jt 


'0 JT 

We conclude by Lemma 6 . 


□ 


5. Proof of Theorem 1 


For the sake of brevity, we are going to prove the case B defined in (7). 
The remaining case B given by (9) requires only minor modifications. We 
only consider the case x > r, being the other case much easier. 

5.1. Approximate problems. Let T > 0 be a fixed parameter. For e G 
( 0 , 1 ) we define the following approximate problem 

(37) dtUe = -eh^-'^^Ue - + xdx{ueB{ue)) + /(Ue), 

where B is given by (7) (or by (9)) and / is given by ( 8 ). The function u 
verifies the initial condition 

Ue(0, x) = Je* uo(x) > 0; 

where 77^* is the heat kernel at time e. Along the lines of Lemma 3 we have 
local existence of solution. Let us define 


and 


Ue{xt,t) = maxue(x,t), 

X 

2M 

ro = 


Ue{xt)' 

We have two possibilities, namely vq > ir or < vr. Assume first that 
ro > TT. Then we have the bound 


(38) 


2M 

— > ||w(f)||L-. 

TT 


On the other hand, if ro < vr, we can apply Lemma 5 with p = 1,71 = A7 
and we get 

A^-'^®n,(xt) > C{M)\\ue{t)\\Y^. 

Using pointwise methods as in [1, 11, 39], we have 
d 


\Ue{t)\\L-° < -eC{Af)\\ue{t)\\lIJ‘+x\\ue{t)\\l^+r\\u, 


dt 

As a consequence, we obtain a uniform-in-time bound 
(39) ||u,(t)||Loc < C(e,x,r,A7). 

Collecting both inequalities (38) and (39), we have the bound 

||w(f)||L°° < C{e,x,r,Jf) + —, 

TT 


.(l-||Uj 
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regardless of the value of rg. 

The continuation criterion implies that ttg is a global-in-time, smooth 
solution (recall that T is arbitrary, but fixed). 

5.2. Uniform estimates. Part 1. We will use in what follows thesis of 
Lemma 4 for in place of u, since it holds for Uf, by the same token as for 
u. 

Let us consider first the case r > x/2. Using this condition, we take 

r 

s = l< -. 

X-r 

Lemma 4 gives the following uniform bounds for all t < T: 

(40) sup \\Ueit)\\Ll <M, 

0<t<T 

(41) / ||ue(s)||^„/2_5,i(is < C(q;,( 5,W, ||t6o||L2,T), 

^ 0 

for arbitrary 0 < J < a/2. 


(42) 


sup ||u£(t)||ii+. < C{r,T,\\uo\\Li+s), 
0<t<T 


(43) 


(44) 


[ \\ue{s)\\l 2 ds < C{r,T, \\uo\\li,x)- 
Jo 

{r{s + l)-xs) [ \\u{s)\\/^‘/ 2 ds <C{r,T,\\uo\\Li+s), 

Jo 


(45) 

for arbitrary 0 < 5 < a/(2 -|- 2s), 

Furthermore, by (33) and following the proof of Lemma 4, we also have 

(46) [ ||ue(s)||^„/ 2 ds < C'(r,r, ||uo||l 2 ,x). 

Jo 

Part 2. Now, let us consider the case r < x/2. Using Lemma 4 and 
taking 

r 


(47) 


s = 


« 1 ), 


X - ?• 

the uniform bounds (40)-(43) hold. We do not have (46), but Lemma 7 and 
(42) imply 

J ll'*^e('®) ll^ya/(2+2s)-«,l+s'^'® — C{a, S, S, ||tt0 ) ^) 

for arbitrary 0<5<a/(2-|- 2s). 
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5.3. Convergence. Due to the uniform bound (43) we have 

dtu,GL\0,T-,H-^) 

e-uniformly bounded. 

Part 1. Now we first consider the case r > x/2. We take 
= H-"^, X = L^, Xq = 


Then, due to the Aubin-Lions Theorem, we have that 

Y = {h: he L\0 ,T-,Xo), dth E L‘^{0,T-, Xi)} 


is compactly embedded in 


L\0,T;X). 


Hence we can extract a subsequence (denoted again by u^) such that 

Ue^uin l2(0,T;H'"/2), 


Ue —)■ u in L^(0, T ; L^). 

This strong compactness is enough to pass to the limit in the weak formula¬ 
tion. As a consequence, n is a global weak solution (in the sense of Definition 
1) of the original system (6), (7) (or (6), (9)) with (8). Its regularity follows 
from l.w.s.c. of norms and uniform bounds (41) - (46). 

Part 2. Now we explain how to adapt the proof for the second case 
r < x/2. Let us first use the assumption a > 1 — s. It implies 


1 1 
2 ^ 1 + s 


1 - - 
2 


consequently, for <5 > 0 sufficiently small, 

y^a/(2+2s)-6,l+s j2 


We take 

Ai = 77“^ A = Ao = vp«/(2+2D-5,i+s^ 


and we apply Aubin-Lions’s Theorem again. The rest of the proof remains 
unchanged. 

Next, let us focus on the assumption a > , where, abusing nota¬ 

tion, we have redefined 

r 

0 < s < -. 

X-r 


It implies 

yYa/(2+2s)-5,l+s . 


We choose 

Ai = 77"^ A = Ao = ^«/(2+2^)-54+s^ 


the Aubin-Lions Theorem implies now that one can extract a subsequence 
(denoted again by u^) such that 

n, ^ u in L2+2"(0 ,r;L^). 

For the limit passage in the most troublesome part —rul we want to use 
additionally (4) of Lemma 4 for Ug, i.e. 

(r(s + 1) - xs) [ Wu.isW+lAs < e'’("+^)*||uo||2t+i- 
Jo 
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For its left-hand-side to be meaningful, we need to choose now in (47) the 
sharp inequality. Consequently, we have 

r 

r 

Ir. Ir^'^ ^ M ^ M. 

and the right-hand-side vanishes as e —)• 0 thanks to strong convergence in 


■ <r\(j)\L'^L^\\Ue-u\\ ^\\Ue + u\\^2j^s+2 

Jo Jt ^ 


s+2 


L^(0,T;L'>+i) and boundedness in L^+^(0, T; 


6. Proof of Theorem 2 

For the sake of brevity, we are going to prove the case of B defined by 
(7). For the other one {B given by (9)), the proof can be easily adapted. 

For a = 2, the proof follows from standard energy methods and Sobolev 
inequality. Hence, let us focus on the case 2 > a. First, let us fix an 
arbitrary T < oo. The proof is based on obtaining the bound 

sup \\u{t)\\L^ < C{a,r,x,uo,T), 

0<t<T 

where the constant C{a,r,x,uo,T) is finite for finite T. 

A computation gives (see [24, 25, 22, 23]) 


(48) 

with 

/(/) := c„ P.V. 


/A“/ = h-{f) + !/(/), 


fcez\{o}' 


\x - yP 


(fix) - f{y)? 

\x — y + 2/c7rP+“ 


dy. 


Let us multiply (6) by u and define v as in (27). We have that 

+ uK^u = ru^{l — u) + X Hv \ . 

(XL \ ^ J 

Now we use identity (48) for the dissipative term. We obtain 

(49) -|- A"(n^) + I{u) = 2rn^(l — u) + x {^u^K^v + dx{u^)A^~^HvJ . 

We intend to control the term 2xtt^A^u with I{u). The rest on the right- 
hand side will turn out to be harmless. We obtain this control in a few 
steps. 

6.1. Auxiliary inequality for primitive function. Let us define 

/ X 

u{z)dz, 

-TT 

where we suppressed the time dependence and 

r 


0 < s < 


X-r 


We have 


\U{x)-U{x-y)\ = 


rx 

/ u{z)dz 

Jx-y 


< lkllL'*+i|yh+^- 
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Hence, using (3) of Lemma 4, we have for any x G T and 0 7 ^ y G T 
(50) < ||n(t)|U.+i|y|^ < e’'iuo||L«+i|y|^. 

\y\ 2 (s+l) 

Notice that the previous equation is a particular case of the splitting 

1 5 

+ 


s + 1 


1+5 1+5/s+l 


We use this particular split 5 = 1, because it reduces the number of pa¬ 
rameters present in the proof and consequently makes our argument more 
traceable. Notice however, that the value of 5 may be in fact arbitrary with 
no consequences to other parameters. This observation will be used in the 
last step to generalize the proof. 

6.2. Towards control of with I{u) in (49). Let us compute for 

a i? < 1 

(tt(x) — u{x — y))^ 


c-^I{u) > P.V. [ 
Jt 


\y\ 


l+a 


-dy 


> 


{u{x) — u{x — y)y 


> 2u 


TT>\y\>R 
2 / 


|y| 


l+a 


dy 




dy — 2u{x) 


• 7T> 


ujx - y) 


dy 


It reads 

(51) 

with 


Ca I{u)>—{R "-vr ")u (x) - 2u(x)(/i + 12 ). 
a 


h := 


u{x - y) 


dy, h ■■= 


u{x - y) 


dy, 


JL>\y\>R ’ J-K>\y\>L 

where L > i? is to be chosen later. Let us examine Ii, l 2 - Using the identity 

u{x - y) = dy{U{x) -U{x- y)), 

we see that 

r U{x,t) -U{x-y) y 

'L>\y\>R 


II — (1 + Ct) 


\y\ 


2+a 


\y\ 


dy 


+ 


U{x,t) -U{x-y) 


\y\l+a 

Consequently by (50) we get 


+ 


U{x,t) -U{x-y) 


R 


|y| 


l-l-a 


-R 


-L 


Ji < (1 + a) 


'L>\y\>R 


e \\uo\\L‘>+yy\^‘'’’+^^ 

|y| 2 +«- 2 Q^ 


dy 


^ ^ ^e^y\uo\\Ls+iL^ 






Let us denote 


D := e 


rt 


^s+lL2G+l). 
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Via D > we arrive at 

rL 


h < (l + a)2 
Hence we can estimate 


R |y| 2 +«- 2 (^ 


dy + 2 


D 


2(s+l) 


+ 2 


D 


^l+“ 2(s+l) 


Co{a,s)D Co{a,s)D 

-^1 — T I S \ 


2(«+l) I 


l+a— 


2(s+l) 


with 


Cq{oi, s) — 2 + 


2(1 +a) 


2 + a — 


2(s+l) 


For I 2 we use Lemma 4 to get 

h < L-(^+“)||u(t)||ii < AA. 

Estimates for /i ,/2 give in (51) 

Ca^I{u) > - 7r“")tt^(x) - 2u(x)(/i + I 2 ) 

a 

> -(i?-“ - 7r-“)u2(x) - 2n(x)L-(i+“)AA 
a 


(52) 


— 2 ( 70 ( 0 ;, s)u(x) 


D 


+ 


D 


2(s+l) ^ 


1+a— 


2(s+l) 


Let us introduce two parameters 7 , m > 0 to be fixed later. Also, let us now 
choose R and L according to 

1/oL 

R-.= ' 

(53) 


Hence 

and 


^70 u{x) + 

D = e''*||uo||L»+iL 2 G+i) := min {e'’*||rto||L«+i, c„} . 


-R " = + "i) > ^rn 


D < Ca- 

Recall that we need L > R. The formula (53) for D fixes the value for L, 
so we need to choose the appropriate value of m such that L > R. 

In the case when min{e^*||rto||s+i, Co} = e''*||tto||s+i, L = 1 and R should 
verify 

f 2 dV/° / 2 


(54) 


R < 


\jam 


< 


70 m 

,rt I 


<1 = L. 


On the other hand, in the case of min{e”||no||s+i, Ca} = Ca, R should verify 


(55) 


2{s+l) 

L=(c„e-’'‘||«„|| 7 „) ■ > 


2 Cg 

'ja m 


Ifcx. 


> R. 


For equations (54) and (55) to hold, we take 

m = — max 11 , (cae“''*||uo||))i+i) 
a7 [ ' 


2q(s + 1) 
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The only parameter which is still free is 7 . Inserting our choice of R (53) 
in inequality (52), we obtain 

2 
a 

+ 


c^^I{u) > ^u^{x) -TT "n^(x) — 2u{x)L 


— 2 u(x)C'o(q;, s) 


D 


I f 2 ^ V 

\ y7Q: u{x)-\-m J 


(l+““ 2 (s+l))/" r^+a- 


2(s+l) 


Using 


L-i < 


2 {£±lij 

e s 


2 (s-H) ’ 


(cq,||moI7s-i-i) 

we arrive at 

I{u) > 7n^(x) — C'i(a)M^(x) — C2{a, uq, ^ “ ^*u(x) 


— tt(x) 


2caCo{a,s)D 

(i+“-27i^)/" +<^3(a,wo,s)e 


l+Q,_5_J. 2(s+1) j 

^ 2 (s+ 1 )T s 


2 D 

70 ti(a;)+m 


where 


Ci(q!) — -vr , 

a 


( 72 ( 0 ;, 110 , 5 ) = 


2 Cn 


(Cq, 11 Roll ^s+l) 


-1 




C3{a,uo,s) = 


2clCo{a,s) 


Notice that for 




1 — a < 


2{s + 1 ) 

the exponent (1 + a — 2 ( 5 + 1 ) )/q =: cr < 3, so we arrive at 

(56) I{u) > 7 it^(x) — C{m, D, a, s, uo)(u^(x) + 1) — c(q!, uq, s)e'^^'’’*’“^*it(x) 
with 

fj < 3. 

6.3. Control of 2xit^A^i; with/(it) in (49). Let us use (56) in (49). Hence, 
we obtain for p = v? 

dtp + A"p + (7 + 2r)pi < 2rp + x (^2pA^v + dxpA^~^Hv^ 

+ C{m,D,a,s,uo){p'^^'^{x) + 1 ) 

(57) +c(a,uo,5)e^(’'’^’“)*V^. 

Since Lemma 1 (part 1) and Definition (16) together with the positivity of 
u, we have 

(58) A^v < u, 
for both operators B{u) considered. 
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As a consequence, equation (57) can be estimated as 
dtp + A"/9 + (7 + 2r — 2 x )/)2 < 2rp + xdxP^^~^Hv 

+ C{m,D,a,s,uo){p'^^‘^{x) + 1 ) 

(59) +c(a,no,s)e^("’^’“)‘VP- 

We are free to choose the value of 7. Let us take 7 = 2x. Now, defining 

p{t) = max p{x, t) = p{xt,t) 

xgT 

we have 

^p{t) + A°‘p{xt) + 2rpi < 2rp + C{m, D, a, s, uo){p''^'^ix) + 1) 
at 

+ c{a, uq, 

Thanks to r > 0 and fi < 3, we get 

C{m, D, a, s, uq){P^^‘^{x) + 1) — 2rp^/‘^ < C{m, D, a, s, uq). 
Therefore 

^p + A"p(xi)+ < 2rp + C{m, D, a, s, uq) + c(a, uq, 
which implies that 

sup ||u(t)||L°° < C{r, s,a,uo,T,m, D) < 00. 

0<t<T 

Notice that T was fixed at the beginning of the proof but it can be taken 
arbitrarily large. 


6.4. Recovering the range of admissible a’s. Observe that the required 
assumptions were 

• s < (for the estimate in Lemma 4) and 

• 1 < 2 (^ ( 56 )). 

In particular, we can take s = r/{x — r) and we get the global existence of 
solution in the range 

(60) a6(l_ir,2(). 


6.5. Generalizing the range of admissible a’s. Observe that the ex¬ 
pression 


1 r 

2x 

in (60) originated in Section 6.1. In particular, we can change (50) as follows 

(61) 

|y| l-f(5 (s+1) 

for an arbitrary 5 > 0. By doing so, we arrive at the range 


(62) 

Since we can take 0 < 5 <C 1, we recover the range as in the assumption. □ 
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7. Proof of Proposition 1 


Define 

Then we have 


dt 


u 


u{t) = maxu{x,t) = u{xt,t). 


+ K'^u{xt) = xu^^v{xt) + ru{l — u). 


Since (58) holds for v defined in (15) and (16), we obtain 


■^u + k°'u{xt) < + ^^(1 “ u). 


Let 


Fo = 


2M 

u{xt)' 


We have two possibilities, namely fq ^ ^ or tq < tt. In the former case we 
have the bound 

2M 

TT 

On the other hand, if fq < vr, we can apply Lemma 5 with p = l,xi = Af 
and get 


k°^u{xt) > 


Ca U 


Thus for ^ < ||u(t)||ioo holds 


2“ AA" 


Cn U 


l+a 


—u < YU^ + mil — u) - -—rr 

dt - ^ ^ ’ 2°‘ 


This is a differential inequality of type 


(63) + BX^{t) - CX^+^{t) 

Let us focus on the case a > 1. We use a blowup argument to obtain 
the global bound. More precisely, let us denote by sq the lower bound 
of values s for which As + Bs^ — < —1. Consider the case when 

X(0) < max(so, ^). Assuming that there exists first time to > 0 such 
that A(0) = max(so, ^), we obtain from (63) that -^X{t) < —1, which 
contradicts the fact that to is the hrst time of equality. As a consequence, 
X(t) < max(so, ^) for all times. The remaining case W(0) > max(so, ^) 
means that as long as X{t) > max(so,^), we can use (63) that gives 
exponential damping ^X{t) < —1. Therefore at certain ti < oo we have 
A(ti) < max(so, ^)- Now we repeat the argument from the previous case. 

Thus, we obtain the global bound 


sup I 

0<i<oo 


In the ’critical’ case a 


u{t)\\L^ < C{a,r,x, ||wo||l°°, IkollLi)- 

= 1, we have for ^ < ||?x(t)||2,cxi 


1 

27rmax(||Mo||LL27r) 

Under the assumption x ~ ^ ~ 27rmax(|Do|| i Ztt) ^ blowup argument 

implies again the thesis. □ 


^ + ru. 


d _ ( 

< \X-r- 
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Appendix A. Pointwise estimates for the fractional laplacian 


Lemma 5. Let h G C'^(T) be a positive function and write h{x*) = max^, h{x) 
I £, 00 . Then, if h verifies the bounds 


£ 00/2 > (h), ||h||£p(Tr) < 7 p, with any 'jp < 




then 




r(l + a) cos((l — a)7r/2) 1 ||/i||^ 


2 

l+ap 


TT 2P“ 

Proof. Let /? > 0 be a constant to be hxed later. We define 

= {7 € B{0,p) I h{x*) — h{x* — rj) > h{x*)/2}, 

and U 2 = B{0,p) — Ui. Notice that if the function is sharp enough, i.e. if 
h{x*)/2 > minj,g5(o,p) then Ui 7^ 0. Starting from our assumption, we 
have 

\h{x*)\P 


so, via positivity 
(64) 


f \h{x* — r])\^ dp > f \h{x*—p)\^dp 

Jr Ju2 


> 


2P 




^7p 

h{x* 


Recalling (14) together with 


^7p 

\h{x*)\ 


< -\U2\. 


c{a) = 


r(l + a) cos((l — a)7r/2) 


TT 


we have 


A"/i(x*) = c(a) Vp.V. [ 

, Jt 


h{x*) — h{x* — 7 ) 


> c(a)P.V. 


T |7 + /c27r|i+" 
h{x*) — h{x* — 7 ) 


dp 


tUi 


\v\ 


l+a 


dp 


. . h(x*) , . 


> c{a) 


h{x* 


2p- 


Let us fix now 


2pi+“ 

( \ 


h{x* 
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thus 


> c(a) 



c{a) h{x*)^+°P 


Finally notice that due to the boundedness of the domain we have to impose 
the restriction 


r = 



< TT i.e. 7p < 


^||/i||loo 


2 


□ 


Appendix B. Sobolev-type inequalities for the fractional 

LAPLACIAN 

We have the following Sobolev-type bounds for the entropy 

Lemma 6. Let 0 < u E LP{T'^) be a given function and 0 < a < 2, 
0 < 5 < a/2 two fixed constants. Then 

• for p = 1 


A"tt(x) \og{u{x))dx, 

• for p = oo 

\\u\\%c /2 < C{(^,d)\\u\\L<=° / A"u(x) log(u(x))iix. 

For the proof of the case p = 1, the interested reader can see [14] by the 
authors &: Luli. The case p = oo was proved by Bae &: Granero-Belinchon 
in [3]. 

In the next lemma we obtain fractional Sobolev bounds in terms of f n^A": 

Lemma 7. Let 0 < u G s <1, be a given function and 0 < a < 2, 

0 < 5 < a/(2 -|- 2s) two fixed constants. Then, 


Iwck/2 —<5,1 — / 

JTd 


0 < 


/ 


A^u{x)u%x)d 


X. 


Moreover, 


\u 


|2+2s 

I ^q:/(2-1-2s) —5,1-l-s 


< (7(0!, s, (5)||u||][|'+s [ A'^u{x)u^{x)dx. 

Jt 


h = [ u^{x)A^u{x)dx. 

Jt 


Proof. Let us define 
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Using (14) and changing variables, we compute 

u{x) — u{r]) 


h = Co 


/ E /“‘W 




X — rj + 2A:7r|^+" 


drfdx 


+ cj RV, 


= C, 


/T JT 


L.^J 




+ C( 


\x - r]\^+°‘ 
u{ri) - u{x) 

u (ri)- - — - , ,; , dridx 

\r} — x + 2A:7rP+“ 

u{ri) — u{x) 


f RV. 

Jt Jt |x-r/|i+“ 


Hence 


h = 


! {u%x)-u^{ri)) 
-'T 


X — r/|i+" 

u{x) — u{rj) 


k^Tj^k^O 


> 0 . 

In particular, using u > 0 
c, 


X — rj + 2A;7rP+" 

u{x) — u{r]) 


drjdx 


f f St \ sf \\U{x) - u{r]) 

■ (^) -^ (d)) |^_^|i+a dvdx 


h > 


Jt Jt 



T Jt Jo 


X - 7?|^+" 

^ c? ,,, , , ,, , u(x) — u(ri) , , 

((\u(x) + (1 — X)u{r])) ) -j—- 11 dr]dx 


dX'^''""''' ' ' ^ 'V“v//y / |^_.^|l+a 

(ar.) =2^ f f ^ _f_ {u{x) - u{r]))^ 

^ ’ 2 AAyo (A«(x) + (1-A)n(r?))i- \x - 

Let us define /? and I 


dXdrjdx. 


I = llul 


Then we compute 

|m(x) — u(r/)|^+^ 


u = 

-II 


a 


-s, 


2 +2s 

u{x) — tt(77)|^+^ 
jj; _ 7y|l + (l + s)/3 


dxdrj. 


I = 


III 

JtJtJo \x - r]\^+i^+^)d 


IT Jt 


where 


\u{x) — u{r])\ |n(x) — u(7?)|^ |Au(x) + (1 — A)n(r/)|(^ 

Q |a; _ ^|0.5+(l+.)/3-a/2 1^ _ ^\0.5+a/2 _ X)u{r,)\i^-s)/2 

m F{x, r], X)G{x, T], X)dXdxdr], 

_ 

u{x) — u{'q)\ 1 


F = 


\x — r/|0-5+a/2 |Att(x) + (1 — A)u(r/)|(i“^)/2 ' 


G = 


|u(x) — u{r])\‘ 


-|Am(x) + (1 — A)u(r7)|^^ 


1^ _^|0.5+(lH-s)/3—a/2 

Consequently 

(66) I < ||T||j;^ 2 (TxTx[ 0,1 ])I|G'IIl 2 (TxTx[ 0,1])- 


dXdxdr] 
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We have via (65) 

II-^IIl2(TxTx[0,1]) “ 
Since 


III 


^ {u{x) — u{r])Y 


dXdxdrj 


< —h. 


\x —-qY+a (Att(x) + (1 — A)u(r/))i“® Cq,s 


^2 \u{x) — |> / N / m1-Q 

^ = 1^ _ ^|l+2(l+s)/3-a 

/ ^ + \u{r])Y+^ ^ |■u(x)|l+^ + |u(r/)|i+* 


_ 7y| l“t"2(l+s)/5 —a 


|a;_y^|l-2(l + s)5 ’ 


thus 


II^IIl2(TxTx[0,1]) — '^)II'“IIli+s 

Estimates for F and G in (66) give thesis. 


□ 
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